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1  Introduction 

Accurate  prediction  of  store  separation  from  internal  aircraft  bays  is  of  great  importance  in 
commercial  and  military  aircraft  design.  The  flow-field  inside  and  over  the  bay  is  very  complex, 
especially  for  supersonic  flows  involving  shock-boundary  layer  interaction,  compressible  shear 
layer,  and  strong  acoustics.  Understanding  the  physics  of  the  flow  in  and  surrounding  the  bay  is 
of  great  importance  for  safe  release  of  stores  and  controlling  the  acoustic  environment. 

The  supersonic  flow-field  within  and  surrounding  the  bay  is  unsteady,  and  accurate  prediction  of 
the  shear  layer  requires  dynamic  simulations  with  advanced  numerical  and  physical  models.  The 
objective  of  this  paper  is  to  demonstrate  accurate  prediction  of  the  effect  of  the  flow  and  a  control 
method  by  employing  a  turbulence  model  that  is  not  cost-prohibitive  and  is  accurate  enough  to 
simulate  such  complex  flow  phenomena.  The  turbulence  model  employed  for  this  study  is  a  two- 
equation  hybrid  Reynolds-Averaged  Navier-Stokes  (RANS)  [1]  and  Large  Eddy  Simulation 
(LES)  [2]  formulation  known  as  the  Large  Eddy  Stress  balanced  (LESb)  model  [3]. 

2  Turbulence  Model 

The  two-equation  RANS  models  use  k  to  represent  all  of  the  unsteady  fluctuations,  while  the  LES 
models  use  k  to  represent  only  the  spatial  average  of  the  fluctuations  within  a  filter  width.  The 
length  scale  Ib  in  the  balanced  LES  model  (LESb)  is  defined  as: 

(2.1) 

and  the  turbulent  viscosity  is  represented  as: 

lu,=pCJs4k  (2.2) 

Eor  unresolved  length  scales  (A  »  h)  the  model  reverts  to  the  standard  two-equation  model. 

However,  for  resolved  length  scales  (A«l^ ),  the  model  reverts  to  the  LES  model,  and  as  the 

resolved  length  scale  goes  to  zero  (full  resolution  of  all  pertinent  scales)  the  model  approaches  a 
Direct  Numerical  Simulation  (DNS)  of  turbulence. 

Implementation  of  the  LESb  model  in  an  existing  RANS  code  with  a  two-equation  model  (and 
the  ability  to  run  time-accurate  simulations)  is  straightforward.  The  existing  model  solves  for  O) , 
which  implies  an  effective  length  scale  of  the  turbulence  and  enters  the  k  equation  as  a  dissipation 
rate  of  the  turbulent  kinetic  energy  (TKE).  To  implement  the  model,  one  should  compare  the 
computed  effective  length  scale  with  the  grid  scale,  and  limit  the  dissipation  term  accordingly. 
Then  one  can  implement  the  filter  as  follows: 

(  ^ 

<y„=max  co, -  (2.3) 

CCA 

The  effect  of  the  limiter  is  to  increase  the  dissipation  of  turbulent  kinetic  energy.  Thus 
the  TKE  is  reduced  from  that  predicted  by  the  traditional  turbulence  model  such  that 
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length  scales  that  are  resolved  do  not  contribute  to  the  Reynolds  stress  terms. 

For  isotropic  grids,  A  is  defined  as 

A  =  {Volume)"^  =  (2.4) , 

where  J  is  the  Jacobian  of  the  grid  transformation.  For  stretched  grids,  the  model  assumes  that 
the  smallest  resolved  eddies  should  be  roughly  isotropic  and  so  must  be  resolved  in  all  three 
coordinate  directions,  and  in  time.  Thus  the  filter  is  defined  as: 

A  =  max((ix,  dy,dz,u^  dt,  4k  *  dt)  (2.5) 

The  limiting  length  scales  based  on  time-step  represent  the  scales  based  on  convection  velocity 
and  SGS  turbulence  respectively.  These  scales  are  included  to  ensure  that  there  is  sufficient  time 
resolution  to  resolve  the  captured  physical  phenomena.  That  is,  if  the  time  steps  are  too  large,  the 
unsteady  phenomena  cannot  be  resolved,  and  the  RANS  model  should  be  used. 

Near  the  wall  in  typical  boundary  layer  grids,  dy  will  be  small  to  resolve  the  mean  shear,  but  dx  is 
typically  very  large.  Thus  the  LESb  model  should  revert  to  the  underlying  two-equation  model 
near  walls.  However,  as  the  boundary  layer  is  traversed  into  grids  that  are  more  isotropic  and 
resolve  the  large-scale  structures,  the  model  should  smoothly  transition  to  an  LES  model. 

3  Grid  and  Solutions 

To  understand  the  impact  of  the  shear  layer  from  an  internal  bay,  we  will  begin  our  study  with  a 
supersonic  flow  over  an  empty  cavity;  then  we  will  investigate  how  to  actively  control  the  shear 
layer  and  suppress  resonance  in  the  bay  by  blowing  with  a  sonic  slot-jet  upstream  of  the  bay.  A 
similar  study  was  carried  out  experimentally  by  Zhuang,  et.al  [4]  and  a  numerical  investigation  of 
approximately  similar  geometry  using  LES  turbulence  model  was  carried  out  by  Rizzetta,  et.al 
[5]. 

A  3D  grid  was  generated  for  a  cavity  (L=12,  D=2.4,  W=2.1  centimeters)  with  an  aspect  ratio  of 
L/D=5.0  The  geometry  contains  the  upstream  plane  to  capture  the  boundary  layer  accurately  and 
side  panels  to  include  the  boundary  layer  interaction  with  the  shear  layer.  The  computational  grid 
contains  2.7  million  points.  To  avoid  a  small  time-step  in  advancing  the  governing  equation  in 
time,  the  grid  spacing  at  the  wall  is  set  to  y^=20.  The  wall  function  [6]  is  invoked  point  by  point 
when  the  y"^  is  greater  than  15.  The  grid  is  shown  for  every  other  point  in  Eigure  1.  The  free- 
stream  flow  conditions  are  Mach=2.0,  Po=2. 17x10^  Pa,  and  To=336  K,  Re/cm=2. 33x10^.  The 
unsteady  cavity  solutions  were  obtained  with  symmetric  and  full  configurations  with  and  without 
control  to  investigate  the  asymmetric  effect  of  the  turbulence  model.  The  plenum  total  pressure 
and  temperature  for  the  slot-jet  are  217.00  psi  and  621.6  ^R. 

The  solutions  were  obtained  using  the  BCED  code  [7].  The  BCED  code  is  an  implicit  general- 
purpose  Euler  and  Navier-Stokes  solver  with  hybrid  unstructured  and  structured/unstructured 
grids.  Numerous  special  algorithms  from  first-order  to  fifth-order  accurate  are  available  for 
structured  grids  applications. 

The  Reynolds  average  Navier-Stokes  equations  with  the  hybrid  turbulence  model  of  Reference  3 
are  used  to  advance  the  governing  equations  in  time.  The  Roe  scheme  with  total  variation 
diminishing,  which  is  second-order  accurate  in  the  physical  domain,  was  employed  in  this  study. 
Since  the  computational  domain  was  split  into  26  zones,  a  global  Newton  with  five  sub-iterations 
was  employed  to  unify  the  zonal  solution  in  time.  A  At  of  1.0  microsecond  was  employed  to 
I  advance  the  solution  in  time,  which  will  provide  one  free^stream  convective  time  scale  in  200 
time-steps. 
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4  Results 

Two  unsteady  solutions  with  the  hybrid  turbulence  model  were  obtained.  The  time-mean  of 
pressure  contours  at  the  center  plane  of  the  cavity  is  shown  in  Figure  2.  The  pressure  in  and  over 
the  cavity  clearly  has  been  suppressed  by  the  slot-jet.  There  is  a  region  of  high  pressure  at  the  aft 
wall  followed  by  a  low  pressure  region  in  the  uncontrolled  case.  The  pressure  in  the  same  region 
for  the  controlled  case  is  lower.  The  pressure  along  the  cavity  wall  at  the  centerline  is  plotted  for 
both  cases  in  Figure  3.  The  pressure  is  lower  at  all  locations  for  the  controlled  case.  However, 
the  lowest  pressure  is  at  x=12.0  inches.  The  time-mean  of  the  stream-wise  velocity  contours  is 
shown  in  Figure  4.  There  is  a  separation  region  upstream  of  the  slot-jet  which  causes  a  secondary 
shock  that  combines  with  the  slot-jet  shock  and  forms  a  lambda  shock.  A  thickening  of  the 
boundary  layer  can  be  observed  upstream  of  the  lambda  shock  as  a  result  of  this  separation. 

Instantaneous  snapshots  of  the  density  gradient  without  and  with  control  at  t=ti  and  t2  are  shown 
in  Figures  5  and  6.  The  shock  structure  due  to  mass  injection  substantially  changes  the  flow 
structure  in  the  shear  layer  and  cavity.  The  shear  layer  is  lifted  off  the  bay,  and  as  a  result  there  is 
reduced  interaction  between  the  shear  layer  and  the  aft  end  of  the  bay.  This  would  substantially 
diminish  the  sonic  fatigue  by  damping  the  feedback  loop  generated  by  the  shear  layer  interacting 
with  the  aft  face  of  the  cavity.  In  the  controlled  case,  the  cavity  leading  edge  shock  is  a  pulsating 
shock.  This  shock  is  formed  due  to  the  entrainment  of  the  flow  across  the  leading  edge  of  the 
bay.  The  iso-surfaces  of  vorticity  magnitudes  are  shown  with  the  density  gradient  in  Figure  7 
with  and  without  mass  injection  at  time  t2.  The  iso-surface  for  the  uncontrolled  flow  shows 
substantially  more  three-dimensionality  in  the  shear  layer  than  the  controlled  flow  shear  layer. 
Figure  8  shows  the  density  gradient  at  the  center  plane  with  several  axial  cuts  along  the  cavity  for 
both  cases  (with  and  without  slot-jet).  The  interaction  of  the  cavity  side  on  the  shear  layer  flow 
and  the  extent  that  the  free  shear  layer  extends  beyond  the  edge  into  the  bay  can  be  seen.  In  the 
downstream  of  the  cavity,  flow  is  being  sucked  into  the  cavity  from  the  side  wall,  while  in  the 
upstream  region,  flow  is  spilling  out  of  the  cavity  over  the  side  wall,  creating  a  weak  vortex. 

The  snapshot  of  the  Mach,  velocities,  and  span-wise  vorticities  are  shown  at  two  time  intervals  in 
Figures  9-12.  The  dynamic  of  the  cavity  flow  is  formed  by  the  large  structures  formed  in  the 
upstream  section  of  the  cavity  that  convects  downstream  as  the  unsteady  shear  layer  oscillates 
over  the  cavity. 

The  pressure  spectra  at  three  different  locations  at  mid-span  on  the  upstream,  center  and  aft  cavity 
walls  are  shown  in  Figure  13.  The  coordinate  locations  are  x=20.3,  25.4,  32.5;  y=-0.3,  -2.4,  -0.3; 
z=0.0  centimeters  accordingly.  The  primary  cavity  tones  are  seen  in  the  1000  to  3000  Hz 
frequencies,  with  primary  peak  at  3000  Hz.  There  are  no  significant  differences  in  the  magnitude 
of  the  frequencies  in  the  upstream  and  cavity  floor.  However,  the  magnitude  of  the  SPL  is  greater 
by  about  20  dB  on  the  downstream  wall  just  below  the  cavity  lip.  The  comparison  of  the  SPL 
between  no  flow  control  and  flow-control  shows  that  the  primary  peaks  have  been  eliminated  and 
the  broadband  spectrum  has  been  reduced  by  about  10  to  20  dB  when  the  control  is  on.  The 
experimental  data  of  Reference  2  shows  about  the  same  order  of  magnitude  reduction  in  the  SPL 
as  we  have  observed  in  the  CFD  study.  However,  in  their  cases  the  pressure  peaks  are  about  20 
dB  below  what  CFD  has  predicted. 

5  Conclusions 

The  hybrid  LESb  turbulence  model  is  a  cost-efficient  approach  to  predicting  complex  high  shear 
flow-field  and  massively  separated  flows.  The  results  indicate  that  the  slot-jet  has  significant 
impact  (20  dB)  in  suppressing  the  acoustic  load  and  the  behavior  of  the  shear  layer. 
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Figure  2  -  Pressure  contours  with  and  without  control  at  mid-span 


Figure  3  -  Pnns  along  the  cavity  floor 


Figure  1  -  Center  plane  grid  at  every  other  point 


to 


Figure  4  -  Time  mean  of  the  stream-wise  velocities  at  mid-span 


5 


3'^  International  Symposium  on  Integrating  CFD  and  Experiments  in  Aerodynamics 
20-21  June  2007 

U.S.  Air  Force  Academy,  CO,  USA 


Figure  5  -  Instantaneous  density  gradient  without  slot-jet  at  time=ti,  t2 


Figure  6  -  Instantaneous  density  gradient  with  slot-jet  on  at  time=ti,  t2 


Figure  7  -  Instantaneous  density  gradient  with  iso-surfaces  of  vorticity  magnitude  at  time=  t2 


Figure  8  -  Uncorrelated  instantaneous  density  gradient  at  mid-span  and  several  axial  flow  planes 
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Figure  9  -  Instantaneous  Mach  contours  without  and  with  AFC  at  time  is 


-0.20 

“-0.80 


Figure  10  -  Instantaneous  stream- wise  velocity  without  and  with  AFC  at  time 


Figure  1 1  -  Instantaneous  normal  velocity  without  and  with  AFC  at  time  is 


Figure  12  -  Instantaneous  span- wise  vorticities  without  and  with  AFC  at  time  ts 
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Figure  13  -Pressure  spectra  at  mid-span  on  the  upstream,  bottom  center,  and  aft  cavity  wall. 
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Problem  Description 


Objective:  Develop  a 
robust  approach  for 
safe  separation  of  store 
from  internal  bay  at 
supersonic  speed. 

-  Problem  description: 
Supersonic  flow  over  a 
cavity 

•  The  shear  layer  over  the 
bay  and  its  interaction  with 
the  downstream  wall 
prevents  the  safe 
separation 


Solution  Approach 


•  Modify  the  supersonic  shear  iayer  by  an  active  or 
passive  approach 

-  An  experiment  with  slot-jet/micro-jet  have  demonstrated 
successful  separation 

-  Unsteady  CFD  solutions  were  obtained  for  empty  cavity  with  and 
without  slot-jet  blowing  to  understand  the  behavior  of  the  shear 
layer 

-  The  cavity  flow  recirculation  changes  due  to  the  blowing 

-  The  shear  layer  is  lifted  due  to  the  blowing 


Re/cm=2. 33x10^ 


Grid  and  Boundaries 


Jet  conditions:  Viscous  wall 

Po=217psi,  1^=621.6  R 


•  A  3D  hyperbolic  structured  grid  generated  to  avoid  zonal 
boundaries  at  the  leading  and  trailing  edge  of  the  cavity. 

•  Grid  stretching:  leading  and  trailing  edge  of  the  bay,  slot- 
jet,  shear  layer  region,  inflow,  and  solid  boundary. 

•  Wall  function,  y^=20 

•  Computational  Grid  =2.7  million 

•  Every  5th  point  shown  for  clarity. 


CFD  Solution  Procedures 


•  Algorithm:  Roe  second-order  accurate  in 
Physical  domain. 

•  Internal  boundaries:  second-order  Roe. 

•  Five  Newton  subiteration  was  performed. 

•  Time  scale:  -t=1 .0x1  O'®  sec. 

•  Prior  to  time-accurate  solution  locai  time-step 
solutions  were  obtained. 

•  Data  collections: 

-  Flow  four  times  crossed  the  cavity  prior  to  data 
collection. 

-  Data  were  stored  at  every  25  microseconds 


Turbulence  Model 

LESb(AIAA-2001-2561) 


LES: 


-  turbulent  viscosity  is  proportional  to  shear  and  filter  width 
squared  (Smagorinsky) 

ju,=piC,Af4^ 

-  turbulent  viscosity  is  proportional  to  square  root  of  the  kinetic 
energy  of  unresolved  scales,  and  the  filter  width  (Kim  &  Menon) 


•  RANS: 


M,  = 


-  Turbulent  viscosity  is  related  to  turbulent  kinetic  energy 

l^t  -  P^Jw  ^  k  =  ^u\  u\ 

-  Define  an  auxiliary  equation  to  obtain  the  length  scale 

•  turbulent  dissipation  rate  (co)  ^1/2 


Balanced  Stress  Model  (LE 
Employed  in  this  Study  (aiaa-2ooi-2561) 


Define  k  as  the  kinetic  energy  of  unresolved  scales 

Use  the  well  defined  k  equation  to  represent  turbulent  viscosity 

m,  =  pCJb^ 

Limit  iength  scaie  based  on  local  grid  resolution 

Ig  =  min(4 ,Cg^)  A » RANS  A « LES 
Use  well  calibrated  2-equation  models  for  high  shear  regions 
Directly  compute  large  scale  unsteadiness  (limit  models  to  small  scales) 

Use  iength  scaie  to  increase  the  dissipation  of  turbuient  kinetic  energy 


r 


co^  =  max 


.1/2 


A 


CO, 


V 


CgCg^J 


Ml 


Ml 


1  = 


03 


CgCO 


or  CO  = 


C  I 

JU  CO 


For  resolved  scales,  k  is  dissipated  -  limiting  turbulent  viscosity 
For  unresolved  scales,  the  2-equation  model  dissipates  high  shear 


Space-Time  Filter 


LESb(AIAA-2001-2561) 


•  Most  LES  formulations  rely  on  a  spatial  filter,  and  it  is  natural  to  set 
the  filter  width  based  on  grid  spacing 

A  =  max(  Jx,  dy,  dz) 

•  RANS  and  URANS  refers  to  a  time  average.  Thus  it  is  natural  to 
think  in  terms  of  what  time  scales  are  being  resolved. 

•  To  model  high  shear,  we  may  need  to  introduce  implicit  operators 
that  allow  a  large  time  step. 

•  we  introduced  time  resolution  into  the  filter  width  to  ensure  that  the 
unsteadiness  is  resolved  in  time  and  space. 

-  Convection  velocity 

-  Turbulent  fluctuation  time  scale 

-  In  the  limit  of  infinite  time  step,  RANS  is  the  appropriate  approximation. 

A  =  max((ix,  dy^  dz^  u  *  dt^  4k  *  dt) 


BCFD  Code 


•  General  purpose  Euler  and  Navier-Stokes  solver 

•  Hybrid  unstructured  and  Structured/  unstructured 

•  Implicit,  multi-zone,  and  overset  grid  capabilities 

•  Numerous  spatial  operators  (1-5*^  for  structured  and  1-2"^  for 
unstructured  grids) 

•  Several  turbulence  models  (S-A,  SST,  k-M,hybrid  RANS/LES 
(LESb,  DES),  Reynolds  stress) 

•  Generalized  chemical  kinetics 

•  Dynamic  memory  allocation,  parallel  (PVM  &  MPI),  platform 
portable,  CFF  and  CGNS  format 


Cavity  Flow  without  Control 

Density  Gradient 


Cavity  Flow  Without  Control 


Density  Gradient 
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Cavity  Flow  with  Control 


Density  Gradient 


Cavity  Flow  Without  Control 


Iso-surfaces  of  vorticity  magnitude  colored  by  Mach  Number 
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Cavity  Flow  With  Control 


Iso-surfaces  of  vorticity  magnitude  colored  by  Mach  Number 


Cavity  Flow  without  Control 


Density  Gradient  Animation  at  Center  Plane 


T=  0.003954 


Cavity  With  Control 


Density  Gradient  Animation  at  Center  Plane 


Cavity  Flow  without  Control 


Vorticity  magnitude  Colored  by  Mach  Number 


Cavity  Flow  With  Control 


Magnitude  of  Vorticity  Iso-surfaces  Colored  by  Mach  Number 


Mean  Pressure  at  Center  Plane 


Mean  Axial  Velocity  at  Center  Plane 
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Pressure  Spectra  at  Mid-span  Along 

the  Centerline 
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Conclusions 


•  Flow  over  an  internal  bay  at  supersonic  flow  with  and 
without  control  has  been  analyzed 

•  The  shear  layer  over  the  bay  without  control  has  a  strong 
interaction  with  the  downstream  bulkhead 

•  The  shear  layer  is  lifted  from  the  bay  due  to  the  slot-jet 
blowing 

•  Developed  a  reliable  and  affordable  numerical  approach 
for  solving  supersonic  flow  over  an  internal  bay 

•  It  is  essential  to  investigate  the  effects  of  passive 
approach  in  controlling  the  shear  layer 

•  Demonstrate  dynamic  separation 


